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ABSTRACT 

We present results of simulations of the distribution of 1809 keV radiation 
from the decay of 26 A1 in the Galaxy. Recent observations of this emission line 
using the Gamma Ray Imaging Spectrometer (GRIS) have indicated that the 
bulk of the 26 A1 must have a velocity of ~ 500 km s 1 (Naya et al. 1996). We 
have previously shown that a velocity this large could be maintained over the 
10 6 year lifetime of the 26 Al if it is trapped in dust grains that are reaccelerated 
periodically in the ISM (Sturner & Naya 1999). Here we investigate whether 
a dust grain velocity of ~ 500 km s -1 will produce a distribution of 1809 keV 
emission in latitude that is consistent with the narrow distribution seen by 
COMPTEL. We find that dust grain velocities in the range 275 - 1000 km s" 1 are 
able to reproduce the COMPTEL 1809 keV emission maps reconstructed using 
the Richardson-Lucy and Maximum Entropy image reconstruction methods 
while the emission map reconstructed using the Multiresolution Regularized 
Expectation Maximization algorithm is not well fit by any of our models. The 
26 Al production rate that is needed to reproduce the observed 1809 keV intensity 
yields in a Galactic mass of 26 Al of ~ 1.5 - 2 M 0 which is in good agreement 
with both other observations and theoretical production rates. 

Subject headings: gamma rays: observations — gamma rays: theory nuclear 
reactions, nucleosynthesis, abundances — supernovae: general supernova 


remnants 
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1. Introduction 

In this paper we present a model to explain the spatial distribution of the Galactic 
emission at 1809 keV. The prediction (Ramaty& Lingenfelter 1977, Arnett 1977) and 
subsequent detection (Mahoney et al. 1984; Share et al. 1985)) of this nuclear decay 
line from the decay of 26 Al has been one of the major successes of theoretical nuclear 
astrophysics and gamma-ray astronomy. The main production mechanism for 26 Al is proton 
capture on 25 Mg and thus production sites must be proton and/or Mg rich. Proposed sites 
include core-collapse supernovae (SNe), Wolf-Rayet stars, asymptotic giant branch (AGB) 
stars, and classical novae (Clayton & Leising 1987; Prantzos & Diehl 1996; Diehl Timmes 
1997). 

The COMPTEL gamma-ray telescope on the Compton Gamma-Ray Observatory 
( CGRO ) has added greatly to our knowledge of this emission. It is an imaging Compton 
telescope with an angular resolution of 3.8° (FWHM) and an energy resolution of roughly 
140 keV (FWHM) at 1809 keV (see e.g. Knodlseder 1997). Maps of the Galactic distribution 
of this radiation show that the emission is strongly concentrated to the Galactic plane 
(Diehl et al. 1997, 1995; Knodlseder et al. 1999a; Oberlack 1997) but the exact extent of 
the emission in Galactic latitude is not well constrained because of the the coarse angular 
resolution of the COMPTEL data (Knodlseder 1997). In order to better understand 
the source (s) of the 26 Al, the overall distribution of this emission has been compared 
with both those expected from simple source models such as exponential disks (Diehl et 
al. 1997; Knodlseder 1997) and emission maps at other wavelengths (Knodlseder et al. 
1999b). Knodlseder et al. (1999b) found that the 1809 keV map was most similar the 53 
GHz free-free emission map obtained using the COBE DMR (Bennett et al. 1992). This 
microwave emission is a tracer of ionized gass in the interstellar medium (ISM) and thus 
of massive stars in the Galaxy. This agreement between the distributions suggests that 
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sources associated with massive stars such as core-collapse supernovae and Wolf-Rayet stars 
are more likely to be the producers of 26 A1 than either AGB stars and classical novae. 

Studies of the kinematics of the 26 A1 have been made possible through the use of high 
resolution germanium spectrometers. Naya et al. (1996) used the balloon-born Gamma Ray 
Imaging Spectrometer (GRIS) experiment to do such a study. They found that the intrinsic 
width of the 1809 keV line during Galactic Center transits was 5.4 (+1.4, -1.3) keV. The 
COMPTEL data does not have the energy resolution to confirm or contradict the GRIS 
line width results. This large width requires that the decaying 26 Al atoms have velocities 
of ~ 540 km s -1 which must be maintained over the ~ 10 6 year lifetime of the 26 Al. 

This is problematic since the Coulomb energy-loss timescale in the ISM for singly ionized 
atomic 26 Al is only ~500 years for a density of 1 cm -3 (see e.g. Sturner & Naya 1999). 
Models to explain this observation have concentrated on binding the 26 Al into high-speed 
dust grains where the dominant energy loss mechanism in the ISM is thought to be direct 
collisions between grains and atoms (Drain & Salpeter 1979). There is mounting evidence 
that such high-velocity dust grains condense within the expanding ejecta of supernovae 
not long after the explosion. Early indications of this possiblilty came from analysis of the 
infrared (IR) light curves of extragalactic supernovae (Dwek et al. 1983) and analysis of 
isotopic abundances in meteorites (Clayton 1982). Final proof came from the optical and 
IR lightcurves of SN 1987A which show a dimming at optical wavelengths while the IR 
emissions brightened (see e.g. Dwek 1988; Dwek et al. 1992). 

Chen et al. (1997) suggested a model in which the initial velocity of the dust grains 
is a consequence of the bulk motion of the ejecta during the early expansion phase. They 
found that while the energy-loss timescale of the grains was significantly longer than that of 
atomic 26 Al, the dust grains must still be injected into a very low density (n H < 10 -2 cm -3 ) 
environment if they are not to lose most of their energy prior to the decay of the 26 A1. 
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Ellison, Drury, Sc Meyer (1997) proposed that shock accelerating dust grains that contain 
26 A1 in the ISM may be an explanation for the GRIS observation. We have previously 
synthesized these concepts into a model for the production and maintenance of a high 
velocity distribution of 26 A1 (Sturner Sc Naya 1999). This model was able to explain the 
observed line widths given that the 26 A1 is bound into fairly large dust grains (~ 10 5 cm), 
the ISM has a mean density near 0.2 cm -3 , and the ISM is permeated by SNR shocks 
with a mean separation of 100-200 pc. We found that the momentum of smaller dust 
grains (~ 10 -6 cm) could not be maintained for the ^ 10 5 year time periods between shock 
encounters. 

The large velocity suggested by the GRIS measurement also has implications for the 
spatial distribution of the 1809 keV emission. As a rule of thumb, a particle moving at 1 
km s -1 will travel ~1 pc in 10 6 years. Thus at the velocity required to explain the line 
width, the 26 A1 will travel ~ 540 pc in its lifetime (assuming non-diffusive motion). This 
distance should not, however, be interpretted as the scale height of the 26 A1 distribution. 
This would only be the case if all of the 26 A1 were ejected perpendicular to the disk of the 
Galaxy in a pipe-like outflow. Here we examine whether the narrow latitudinal distribution 
of 1809 keV emission on the sky as measured by COMPTEL can be reconciled with the 
distribution expected from high velocity 26 A1 assuming a more realistic radial outflow from 
the sources. We make no attempt to fit the longitudinal profile of the 1809 keV emission 
which has been examined in detail by Knodlseder et al.( 1999b) and any analysis requires 
the use of complicated source models. In Section 2 of this paper we discuss our model for 
calculating the Galactic distribution of 26 A1 and the 1809 keV sky map. In Section 3 we 
compare the results of our model calculations with the COMPTEL results and discuss what 
contribution the upcoming INTEGRAL mission can make to the field. 
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2. Model 

We have produced a simple model to calculate the distribution of the 1809 keV emission 
which has three parts, 26 A1 production, 26 A1 transport, and sky map generation. In this 
model the 26 A1 is generated within the Galaxy with a production function 

HP 

-jy-fafyP 2 *) = /(ft) g{z s), (1) 

where we have assumed that it is a separable function of the cylindrical coordinates p s and 
z s (azimuthal symmetry is assumed). In this work we will assume for simplicity that most 
of the 26 A1 in the Galaxy is produced in Type II SNe. Stecker & Jones (1977) suggested a 
functional form for the SNR surface density 

where p 0 is the distance from the Sun to the Galactic Center which we take to be 8.5 kpc 
(we assume that z 0 = 0 kpc). Case & Bhattacharya (1998) have re-calibrated the E — D 
relationship for Galactic SNRs and have used this result to derive the values a = 2.00 ±0.67 
and /? = 3.53 ± 0.77 which we use here. Using these parameter values, the surface density 
rises from zero at p s = 0 kpc to a peak at p s = 4.8 kpc from which point the density 
decreases (see Figure 1). 

The distribution of SNe perpendicular to the disk is generally assumed to be that of 
an exponential disk 

g(z s ) = exp(-|z s |/z 0 ) • (3) 

The value for the supernova scale height that we use in our calculations is z 0 = 90 pc, 
assuming it is the same as for massive stars ( M ^ 3M 0 ) (Miller & Scalo 1977; Bahcall 
& Soneira 1980; Ferriere 1993). We will treat e Q in equation (1) as a free parameter to 
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adjust the normalization of the emission to fit the observations. We will then calculate the 
Galactic mass of 26 A1 using this normalization. 


Once produced, we assume the 26 A1 flows radially away from the source locations 
(p s ,<t>s,Zs) at constant velocity v 0 . During the radial outflow, the 26 A1 decays with a 10 6 
year time constant, r 0 . The density of 26 A1 at position (p, 0, z) due to production in volume 
dV s about the position (p s , 0 s ,z s ) is given by 

dn(p, 0, z, ft, </>., z,) = — exp (-<//*. r„) (4) 


where d = d(p , 0, 2 , p s , 0 S , z s ) is the distance from the source location to the reference point 
(p,<£,z). The three terms in equation (4) represent the production, decay, and outflow of 
the 26 Al. The total density at position (p, 0, z) due to emission throughout the Galaxy is 
just 

N(p, (j), z) = / dn(p, (j), z, p s , 0s, z s )- ( 5 ) 

^Galaxy 


and the emissivity of 1809 keV photons at point (p, 0, z) can then be written as 
e(p, 0, 2 ) = N (p, 0, z)/r 0 . 


The observed gamma-ray flux from Galactic coordinates ( l,b ) is then just the integral 
along the line-of-sight 


dF 

dA dt dfl 


(‘■»= r ^ r2 dr = h r r) dr ' 


(6) 


where r is the radial coordinate along the line-of-sight and 

p(l,b,r) = \fp 0 + r 2 cos 2 b - 2 rp Q cos b sin l (7a) 

z(b,r) = rsinb. (^) 

There is no attenuation term in equation (6) since the Galaxy is transparent to 1809 keV 
radiation. Thus this model has 7 parameters: e 0 , a, /?, z 0 , Po, 2 G , and v 0 , all of which are 
constrained by observations. 
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We have used the approximation in this model that the dust grain motion is convective, 
not diffusive. We have previously shown that dust grains must interact with shocks in 
the ISM every 100 - 200 pc along their trajectory in order to maintain a velocity of ~ 

500 km s" 1 (Sturner & Naya 1999). The distance a dust grain must travel in the ISM for 
collisions with gas to appreciably change the grain direction (and thus the energy) will 
be similar to this distance if there is to be an equilibrium between shock acceleration and 
collisional energy losses. The mean distance a grain with velocity v a will travel in time t 
is given by L = y/Xv Q t. Using 100 - 200 pc for the mean-free-path A, the diffusion length 
in the 10 3 * * 6 year lifetime of the 26 A1 is 226 - 320 pc for v 0 = 500 km s -1 as compared to 
the 500 pc distance along the particle path. We can then define a mean grain velocity, 

{v) = L(t 0 )/t 0 = \J Xv Q j t 0 = 226 - 320 km s' 1 . This diffusion length and mean velocity 
should be considered a lower limit as the dust grain mean-free-path will grow as they move 
farther from the Galactic plane. In the next section we discuss the results of this model for 
grains with velocities of 100, 275, 500, and 1000 km s 


3. Results 

In our simulations we have set all parameters to those values given in Section 2 except 
for varying the particle velocity from 100 km s" 1 to 1000 km s" 1 to illustrate how the results 
change with the dust grain velocity. We chose to include results for v a = 275 km s ~ ( v ) 
to approximate diffusive motion of the grains with a true space velocity of 500 km s . In 
Figure 2 we show the 1809 keV emissivity as a function of Galacto-centric radius p and 
height z calculated from equation (5). These plots show a broad peak near p = 5 kpc that 
mimics the 26 A1 production function. We find, as expected, that as the grain velocity is 

increased, the emissivity in the plane decreases while it increases at larger z and the ‘‘hole” 

in the production function at the Galactic Center becomes less evident. 
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We used these results to calculate skymaps of the 1809 keV emission in the range 
—90° < l < 90° and -30° < b < 30°. These maps were then convolved with a 3.8° FWHM 
Gaussian to approximate the effect of the COMPTEL point spread function (PSF) on the 
image. In Figure 3 we show normalized contour plots for the convolved and unconvolved 
emission maps for models with v a = 275 and 500 km s -1 to illustrate the degree to which 
our approximate COMPTEL PSF alters the skymap. The effect of the broad PSF on the 
contour plots is significant. We find that much of the emission in the Galactic plane is 
“transport ed : ’ to higher latitudes. We wish compare these skymaps to the COMPTEL data. 

Knodlseder et al. (1999a) compared several algorithms for performing imaging analysis 
of COMPTEL 1809 keV data, including maximum entropy (ME), Richardson-Lucy (RL), 
and multiresolution regularized expectation maximization (MREM). ME reconstructions 
were very sensitive to statistical noise and tended to produce artificial hot spots that 
may be indistinguishable from real point sources. MREM analysis produced a smoother 
map but at the expense of suppressing weak features. They found that MREM yielded a 
latitude profile (averaged over —90° < l < 90°) with a FWHM = 7.5° while the RL and ME 
methods yielded somewhat wider profiles with widths ~ 8°. We have used our convolved 
model skymaps to generate similarly averaged latitude profiles for the emission. In Figure 4 
we compare the model profiles with the MREM, RL, and ME profiles of Knodlseder et al. 
(1999a). This comparison consisted of renormalizing the model profiles so that the peak 
flux for a given model matched the peak flux of a particular reconstructed map. We have 
also shifted the model profiles slightly in b so that they were approximately centered on 
the data. It should be noted that these are not mathematical fits to the data using x 2 or 
any other statistic. They are only qualitative comparisons to the data since our purpose in 
this paper is ascertain whether high-velocity 26 Al is consistent with the observed spatial 
distribution of 1809 keV emission and not to deduce the exact parameters of this somewhat 
simplified model that best fit the data. 
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The three data reconstruction methods produce significantly different profiles. There is 
more than a factor of 2 difference between the peak flux in the MREM profile and the ME 
profile. The fluxes in the wings of the RL amd ME profiles are similar but the MREM profile 
has has significantly less emission at |6| > 5°. We find that the MREM reconstruction is 
poorly fit by all models. This reconstructed map has the brightest emission in the Galactic 
plane while having the dimmest emission at high latitudes. The v a = 1000 km s -1 profile 
fits the data well for —2.5° ^ b ^ 2.5° but it overproduces emission at higher latitudes by a 
large margin. The other models produce less high latitude emission but do not match the 
lower latitude emission. The RL map is best fit by the model with v 0 = 500 km s l . The 
275 km s -1 emission profile is significanty narrower than the core of the COMPTEL profile 
but has the correct emission level at higher latitudes while the 1000 km s -1 emission profile 
fits the COMPTEL data reasonably well in the latitude range -3° <$ b & 3° but the emission 
is too bright at higher latitudes. The ME map is best fit by the 1000 km s -1 model. This 
model fits the data reasonably well in the region —4° ^5^4° with a small excess at higher 
latitudes. The higher lattude emission is best fit by the 500 km s -1 emission. Thus we find 
that while the model with v 0 = 1000 km s _1 can fit all the reconstructed maps at the lowest 
latitudes, the models with v Q = 275 and 500 km s -1 fit the RL and ME reconstructed data 
for \b\ ^ 5° but tend to underpredict the lower latitude emission. Considering the simplicty 
of the model, including mono-energetic dust grains, the results are quite good. Thus we 
find that the spatial distribution of the 1809 keV emssion measured by COMPTEL is not 
in conflict with the GRIS measurement of high velocity 26 Al. 

One way to possibly improve the fit of the model to the data is to use a somewhat larger 
scale height, z Q , for the production function coupled with grain velocity ~ 275 — 500 km s -1 . 
The larger scale height will broaden the lowest latitude emission while the moderate velocity 
(in the context of this paper) will act to keep the higher latitude emission at a lower level. 
A somewhat larger scale height can be justified given that the 90 pc supernova scale that 
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we use is in our calculations is really a lower limit taken from the scale height of massive 
stars. There is also some evidence that massive stars exist (and thus supernovae must 
occur) at significant distances above the plane of our Galaxy. As we have previously stated, 
Knodlseder (1999b) found that the 53 GHz free-free emission map obtained by the COBE 
DMR (Bennett et al. 1992) produced the best fit to the COMPTEL 1809 keV map when 
comparing 31 emission maps with energies from radio through gamma ray. This free-free 
emission is a tracer of ionized gas and thus the O and B stars in the Galaxy. The ionized 
hydrogen gas in our Galaxy can be described using four components, compact and diffuse H 
II regions which are well confined to the disk, as well as two low-density diffuse components 
with significantly larger scale heights (Lockman, Pisano, & Howard 1996; Taylor & Cordes 
1993). The diffuse and compact H II regions have similar vertical distributions with a 
FWHM ^ 50 pc (Lockman et al. 1996). The two diffuse components have vastly different 
distributions. The higher density component (n = 0.1 cm -3 ) dominates the inner Galaxy 
and has a FWHM ~ 260 pc while the lower density (n = 0.019 cm" 3 ) outer Galaxy 
component has a much larger FWHM ~ 1540 pc (Taylor & Cordes 1993). The origin of 
this low density component is poorly understood (Lockman et al. 1996). 

The model normalization calculated for the fits in Figure 4 can be used to calculate 
the 26 Al mass in the Galaxy which we list in Table 1. Those models that best fit the data 
produce ~ 1.5 - 2 M© of 26 Al in the Galaxy given the ~ 1.04 x 10 6 year mean lifetime. 
These masses compare favorably with those calculated from observations and theoretical 
models. Mahoney et al. (1984) calculated the Galactic mass of 26 Al to be ~ 3 Mg. 
Knodlseder (1997) calculated the mass to be 1.6 - 2.2 Mg using exponential disk models. 
Much smaller values for the Galactic mass result if much of the 26 Al is produced locally (£ 
1 kpc of the Sun). Diehl et al. (1995) find that if there are 5-10 local sources each with a 
flux of 3 x 10~ 5 photons cm -2 s -1 , the Galactic mass of 26 Al would be reduced to ~ 1 M 0 . 
Timmes et al. (1995) have calculated the Galactic production rate and mass of 26 Al using 
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models for the production of 26 A1 in Type II SNe as function of progenitor mass as well 
as stellar mass and birth-rate functions. They found that Type II SNe should produce 
2.0 ± 1.0 M g Myr -1 yielding a Galactic mass of 26 A1 of 2.2 ± 1.1 M 0 in agreement with our 
estimate. We do not suggest this is proof for a Type II SN origin for the 26 Al. The required 
2 M© of 26 Al can in priciple be produced by any of the sources mentioned in Section l(see 
e.g. Diehl & Timmes 1997) and the uncertainties in the true 3-dimensional distribution of 
the 26 Al make it difficult to rule out any type of sources. 

The SPI germanium imaging spectrometer on INTEGRAL may help resolve the 
question of the origin of 26 Al in our Galaxy (see e.g. Vedrenne et al. 1999). One of the 
prime objectives of SPI is to map the Galactic Center and the Galactic plane at both 511 
keV and 1809 keV. SPI combines an energy resolution at 1809 keV that should be at least 
as good as the 3.4 keV resolution of GRIS with a 2.0° angular resolution which is superior 
to COMPTEL. The SPI detector array will also have roughly twice the area and volume as 
the GRIS array. Thus the superior quality and quantity of the data should allow for SPI 
to confirm or contradict the GRIS line-width results. The good angular resolution and the 
high sensitivity of SPI to line emission should allow for the determination of line width as 
a function of position on the sky. This used in conjunction with the overall map of the 
emission may allow for a determination of the major source(s) of the Al. 
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Table 1: 

Model Galactic Mass of 

26 A1 


Vo 

Reconstruction 

e 0 

Mass of 26 A1 

(km s" 

-1 ) Method (10' 

- 25 cm" 3 s" 1 ) 

(Mo) 

100 

MREM 

8.95' 

1.82 

100 

RL 

6.25 

1.27 

100 

ME 

4.28 

0.87 

275 

MREM 

11.0 

2.23 

275 

RL 

7.68 

1.56 

275 

ME 

5.26 

1.07 

500 

MREM 

13.6 

2.76 

500 

RL 

9.48 

1.93 

500 

ME 

6.49 

1.32 

1000 

MREM 

18.8 

3.82 

1000 

RL 

13.1 

2.66 

1000 

ME 

8.98 

1.82 
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Figure Captions 

Figure 1: Plot of the SNR surface density finction given in Equation 5. The parameter 

values are those used throughout our models: = 8 - 5 k P c > Q = 2 -°> anci ^ = 3 ’ 53 ’ The 

peak between 4.5 and 5.0 kpc leads to the peak in the emissivity plots in Figure 2 and the 

peaks seen in the emission contours in Figures 3a-d. 

Figure 2: Here we show the 1809 keV emissivity as a function of Galacto-centric radius for 
v 0 = 100, 275, 500, and 1000 km s -1 and z = 0 (solid), 250 (dashed), 500 (dotted), and 1000 
pc (dot-dot-dot-dashed). We have chosen e 0 = 1.0 x 10“ 24 cm" 3 s" 1 (see Table 1). Note 
that as the dust grain velocity increases the emissivity decreases at small z but increases at 
large z, as expected. 

Figure 3: Here we show contour plots of the observed 1809 keV flux for the models with 
v Q = 275 km s -1 and 500 km s' 1 . We show contour plots for maps that have been convolved 
with the 3.8° FWHM PSF (3b and 3d) and those that have not (3a and 3c). The plots for 
each velocity have been normalized by dividing by the peak flux in the unconvolved map. 
This normalization shows the decrease in the observed peak flux that results from having 
a large instrumental PSF. The peaks seen at ~ ±20° result from the peak in the A1 
production function at 4.8 kpc. 

Figure 4: Here we compare longitude-averaged (—90 < l < 90) latitude profiles for the 
multiresolution regularized expectation maximization (MREM) (a). Richardson-Lucy (RL) 
(b), and maximum entropy (ME) (c) reconstructed COMPTEL 1809 keV emission maps 
(solid) with our model with v 0 — 100 km s -1 (dotted), 275 km s 1 (short dashed), 500 km 
s -1 (dot-dashed), and 1000 km s" 1 (long dashed). We have normalized the model profiles 
such that the peak fluxes equal those for the data. We have also slightly offset the flux 
peak to negative values of b to center the model profiles on the data. Note that none of the 
models can reproduce the contrast between the Galactic plane and the higher latitudes seen 
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in the MREM reconstruction. The 1000 km s _1 model best fits the data at lowest latitudes 
while the 275 km s -1 and 500 km s -1 models best fit the emissio for \b\ > 5°. 















